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Abstract 

Analytical solution of a-dynamo equation representing strongly 
torsioned helical dynamo is obtained in the thin twisted Riemannian 
flux tubes approximation. The a factor possesses a fundamental con- 
tribution from torsion which is however weaken in the thin tubes ap- 
proximation. It is shown that assuming that the poloidal component of 
the magnetic field is in principle time- independent, the toroidal mag- 
netic field component grows very fast in time, actually it possesses a 
linear time dependence, while the poloidal component grows under the 
influence of torsion or twist of the flux tube. The toroidal component 
decays spatially with as r~ 2 while vorticity may decay as r -5 (poloidal 
component) where r represents the radial distance from the magnetic 
axis of flux tube. Toroidal component of vorticity decays as r _1 . In 
turbulent dynamos unbounded magnetic fields may decay at least as 
r 3 . PACS numbers:02.40.Hw-Riemann geometries 
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I Introduction 



Despite of the success of the application of the numerical simulations to the 
dynamo problem [1] in plasma astrophysics [2] and in the stretch-twist-fold 
(STF) Vainshtein-Zeldovich [2] mechanism, recently new dynamo analytical 
solutions have been found [3] by using the conformal mapping technique in 
Riemannian manifolds from old Arnold cat dynamo metric [2]. Earlier T. 
Kambe [4] found simultaneous vortical and magnetohydrodynamic (MHD) 
solutions. In this paper an helical dynamo [8] solution of self-induction is 
obtained in vortical strongly torsioned thin twisted magnetic flux tubes in 
Riemannian space [5] where the MHD equations are linear in the magnetic 
field and nonlinear in the velocity flow. Assuming that the poloidal is time- 
independent the toroidal component of the magnetic field grows fast in time, 
actually it grows linear and not exponential. The rate of growing of the 
toroidal component depending on the inverse squared of the radial distance 
of the magnetic axis which possesses curvature and torsion. Recently, Hanasz 
and Lesch [6] have used also a conformal Riemannian metric in £ 3 to the 
galactic dynamo magnetic flux tubes. Pioneering work on the magnetic flux 
tubes as dynamos was done earlier by Schussler [7] , however in his work tubes 
were untwisted and straight. The main advantage of the investigation of the 
isolated flux tube dynamo is that one is able to investigate the curvature 
and twist contributions of the tube to the dynamo action. Twist is actually 
related to the torsion of the magnetic axis of the tube, which makes the words 
strong torsion equivalent to strong twist, which physically is important to 
the twist-kink relation investigated by Alfven [8]. Helical dynamo here is 
understood as the one where the flow describes a circular helix where torsion 
and curvature are constants and equal. The paper is organized as follows: 
In section II the dynamo solution in the Riemann metric representing flux 
rope (twisted tubes) is obtained. In section III the approximate solution is 
presented. In section IV conclusions are given. 

II Helical dynamos in Riemannian space 

In this section we shall be concerned with solving the MHD equations in the 
curved coordinates of a thin twisted magnetic flux tube of Riemann metric 
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which represents a Riemannian line element 

ds 2 = g ij dx i dx j (2) 

if the tube coordinates are (r, 9r, s) [1] where 9(s) = Or — J rds and r is the 
Frenet torsion of the tube axis, K(s) is given by 

K\s) = [1 - rK(s)cos9(s)} 2 (3) 

Since we are considered thin magnetic flux tubes, this expression shall be 
taken as K ~ 1 in future computations. In curvilinear coordinates the Rie- 
mannian Laplacian operator V 2 [1] is 

V 2 = ^[y/ggVdjl (4) 

where dj := and g := detgij where g^ are the covariant components of 
the Riemann metric of flux rope. Let us now start by considering the MHD 
field equations 

V.5 = (5) 

d -> 

—B = Vx[aB]=aVxB + VaxB (6) 

called the a-dynamo equation [?]. Sometimes the a :=< v. Co > parameter 
is constant but here we shall be considering the more general case where 
it depends on the radial and poloidal coordinate. Here uj := Vxt? is the 
vorticity of the dynamo flow. Equation (3) represents the self-induction 

— * — * 

equation. The vectors t and ft along with binormal vector b form the Frenet 
holonomic frame, which obeys the Frenet-Serret equations 

? = ku (7) 

n' = -Kt + rb (8) 
b' = -rn (9) 

where the dash represents the ordinary derivation with respect to coordinate 
s, and k(s, t) is the curvature of the curve, where k = R~ l . Here r represents 
the Frenet torsion. The gradient operator is 

V-f— e-— e— (10) 

ds 60 r 89 6r dr 
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Now we shall consider the analytical solution of the self-induction magnetic 
equation which represents a non-dynamo thin magnetic flux tube. Before 
the derivation of this result is obtained, we would like to point it out that 
it is not trivial, since the Zeldovich antidynamo theorem states that the two 
dimensional magnetic fields do not support dynamo action. Here, as is shown 
bellow, the flux tube axis possesses not only Frenet curvature, but torsion 
as well, and this last one vanishes in planar curves. The magnetic field does 
not possess a radial component and the magnetic field can be split inti its 
toroidal and poloidal components as 

B(r,s,t) = B e (t,r,9(s)) + B a (r)t (11) 

Now let us substitute the definition of the poloidal plus toroidal magnetic 
fields into the self-induction equation, along with expressions 

ee = —nsinO + bcosO (12) 

and 

e r = ncosd + bsinO (13) 

d t ee = uige r — d t nsin6 + d t bcos9 (14) 

Considering the equations for the time derivative of the Frenet frame given 
by the hydrodynamical absolute derivative 

X = d t X + [v.V]X (15) 

— * — * — * 

where X = (t, n, b) is used into the expressions for the total derivative of 
each Frenet frame vectors 

t = d t t+[K'b-KTn] (16) 

ft = Krt (17) 

b = -Kt (18) 

therefore leading to the following values of respective partial derivatives of 
the Frenet frame 

d t t = -tk[1 - nr' 2 ^]n (19) 
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d t n = tk[1 — 




d t b=KT~ 1 ^-n (21) 

— * 

where use has been made of the hypothesis that b = or n'(t, s) = 0, which 
means that the curvature only depends on time. Substitution of these vec- 
torial expressions into expression (14) yields 

d t eg = —uj e e r + 7[r sm^(l — r )i + T cos9n — r sin6b] (22) 

where 7 := (v s — ^t^vq). The other dynamical equation for the Frenet 
holonomic frame is 

dte*e = [rosinOt — [9$ + T ]cos9n — [ujg + To\sin9b] (23) 

note that in the mean field dynamo case, where v = v{B) , equation (6) is 
an eigenvalue problem equation. Dynamo operators and eigenvalue of dy- 
namos in compact Riemannian manifolds have been previously investigated 
by Chicone and Latushkin [9] . Substitution of previous equations into equa- 
tion (6) and splitting these equations along the components of the Frenet 
frame (t, n, b) yields the following three scalar equations 

d t B s + \sin6 1 - r ]r B e = d r (aB e ) (24) 

—d t B e sin9 — B 9 uqcos6 + (cos9B e j — t B s )t = AcosO — CsinO (25) 
d t B e cos6 — B e {u)e + T )sin6 = AsinO + CcosQ (26) 
where functions A and C are given respectively by 

A :— ^-{dea)B s — (d s a)B e — ad s B e (27) 

and 

C := (d r a)B s + -d r B r (28) 
r 

In the next section we shall solve find an approximate solution for strong 
torsioned a — dynamos tubes. Helical dynamo hypothesis of Ko = constant = 
To has been taken throughout these computations. 
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Ill Analytical solution of helical dynamos in 
flux tubes 

From the solenoidal equation for the magnetic field one obtains 

d s B e = B T O 2 rsin6 (29) 

one obtains the value for the poloidal component as 

Bq = B exp(T rcos9) (30) 

with this expression and the expressions for the value of a which can be 
obtained as 

a = vqUq + v s u s (31) 
To compute this important factor a one computes the vorticity co as 

u) r = —d s ve = KoTorvesin9 (32) 

where we have used in this equation the physical assumption of the incom- 
pressibility of the dynamo flow. The remaining vorticity expressions are 

uj e = -d r v s (33) 

u s = d r v e + -v e (34) 
r 

This allows us to obtain the following value for a 

a = -voexp(—rTocos9) (35) 
r 

By making use of thin tube approximation where we consider we are close to 
the magnetic flux tube axis (r = 0), one may use the following approximation 

a = -vq (36) 
r 

which yields 

d s a = ~^v (37) 



6 



which upon substitution in equation (refl) yields 



v 0l 



d t B s = [to 2 - ^]B e (38) 
d t B s = -^B (39) 



which simply solves to 

R. = -I 

„2 



B s = -^B t (40) 



which shows that the toroidal field grows linearly in time which close to the 
magnetic axis may represent a very fast dynamo. Comparison between the 
poloidal and toroidal field one obtains 

B s r A 

In regions not so close to the magnetic axis where torsion may dominate one 
obtains an alternative solution in which however the magnetic toroidal field 
still grows in time as 

B s ^t H (42) 

Substitution of these equations in the remaining dynamo equation yields the 
poloidal component of vorticity as 

U « = ( 43 ) 

which implies that 

and the other vorticity components are 

u r ^i — T 2 rvoexp(rorcos6) (45) 

and u> r = at the magnetic axis. Finally 

u s = —exp(T rcos6 (46) 

which decays slower as we go away from the magnetic axis of flux tube 
dynamo. 
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IV Conclusions 



In conclusion, we have used an approximate method of strong torsion to find 
near the magnetic axis of the dynamo flux tube to analytical (non-numerical) 
solutions of a-dynamo equation. Previously Ruzmaikin et al [10] have found 
solutions of turbulent dynamos where unbounded magnetic field could decay 
at least as r~ 3 which is distinct to the bounded case of twisted magnetic flux 
tube helical dynamo we have found here. The general equations of the helical 
dynamo found in section II could be used to find out more general solutions, 
with less degree of approximation that we used here, for example by letting 
the magnetic poloidal component vary with time and also considered the case 
of untwisted tubes. This can be done elsewhere. 
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